An SU(4) Model of High- Temperature Superconductivity and Antiferromagnetism 



O 

o 
o 

(N 

o 

Q 

\o 

o 
o 

or 

CZ3 
•4-3 



I 

-a 
c 

o 
o 



> 

00 
00 

o 

o 
o 

-t— > 
c3 



I 

-a 
c 

o 
o 



Mike GuidryW, Lian-Ao Wu< 2 \ Yang SunW and Cheng-Li Wu^ A \ 
^Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996 
''^Department of Physics, Jilin University, Changchun, Jilin 130023, PRC 
^Department of Physics, Chung-Yuan Christian University, Chungli, Taiwan 320, ROC 
^National Center for Theoretical Science, Hsinchu, Taiwan 300, ROC 
(February 1, 2008) 

We present an SU(4) model of high-temperature superconductivity having many similarities to 
dynamical symmetries known to play an important role in microscopic nuclear structure physics 
and in elementary particle physics. Analytical solutions in three dynamical symmetry limits of this 
model are found: an 5*0(4) limit associated with antiferromagnetic order; an SU(2) x SO(3) limit 
that may be interpreted as a d-wave pairing condensate; and an SO(5) limit that may be interpreted 
as a doorway state between the antiferromagnetic order and the superconducting order. The model 
suggests a phase diagram in qualitative agreement with that observed in the cuprate superconduc- 
tors. The relationship between the present model and the SO(5) unification of superconductivity 
and antiferromagnetic order proposed by Zhang is discussed. 



There are compelling arguments that the mechanism 
leading to high-temperature superconductivity does not 
correspond to ordinary BCS s-wave pairing. Although 
experimental evidence implicates singlet (hole) pairs as 
the carriers of the supercurrent, the interaction leading 
to the formation of the singlet pairs appears not to be 
the traditional lattice phonon mechanism underlying the 
BCS theory, but rather seems to be a collective electronic 
interaction. Furthermore, the pairing gap is anisotropic, 
with nodes in the k x -k y plane strongly suggestive of d- 
wave hybridization in the 2-particle wavefunctions, and 
the mechanism responsible for superconductivity in the 
cuprates is thought to be closely related to the unusual 
antiferromagnetic (AF) insulator properties of their nor- 
mal states. 

Contrary to the case for BCS superconductors, the for- 
mation of Cooper pairs and the formation of a supercon- 
ducting (SC) condensate of those pairs in high-T c com- 
pounds may be distinct, with pair formation correspond- 
ing to a higher temperature scale than the condensation 
of the pairs into the SC state. That is, there appear to 
be at least two distinct energy scales associated with the 
formation of the high-temperature SC state. This is rem- 
iniscent of grand unified theories in elementary particle 
physics, where qualitatively different physical phases re- 
sult from a hierarchy of symmetry breakings occurring 
on different energy (temperature) scales. This finds its 
most natural explanation in a Lie group structure that is 
broken spontaneously (and perhaps explicitly) down to 
subgroups at different characteristic energy scales. 



I. DYNAMICS AND SYMMETRIES 

Such observations argue strongly for a theory based 
on continuous symmetries of the dynamical system that 
is capable of describing more sophisticated pairing than 
found in the simple BCS picture (which is described by 



a single complex order parameter), and capable of unify- 
ing different collective modes and phases on a equivalent 
footing. Then such fundamentally different physics as an- 
tiferromagnetic order and superconductivity can emerge 
from the same effective Hamiltonian as concentration 
variables (e.g., doping parameters) are varied. 



A. Fermion Dynamical Symmetries 

For approximately the same period of time that the 
high-T c compounds have been known, techniques based 
on dynamical symmetries in fermion degrees of freedom 
that are capable of satisfying the preceding conditions 
have been in development in the field of nuclear struc- 
ture physics. There it has proven fruitful to ask the fol- 
lowing questions: what are the most important collective 
degrees of freedom in the low-lying spectrum of complex 
nuclei, what are the microscopic many-body quantum 
operators that create and annihilate these modes, and 
what is the commutator algebra obeyed by this set of 
operators. 

Systematic investigation of these questions has led to 
strong confirmation of the following set of conjectures 
about the nuclear many-body system: (1) Low-lying col- 
lective modes are in approximate one-to-one correspon- 
dence with dynamical symmetries in the fermion degrees 
of freedom. (A system possesses a dynamical symmetry 
if it has a Hamiltonian that can be expressed as a poly- 
nomial in the Casimir invariants of a subgroup chain.) 
(2) A dynamical symmetry associated with low-lying col- 
lective modes is associated with a Lie algebra and its 
subalgebras that are formed from a set of fermion opera- 
tors closed under commutation. (3) Different dynamical 
symmetry subgroup chains arising from the same high- 
est symmetry group are associated with fundamentally 
different phases of the theory. These dynamical sym- 
metries are characterized by different collective modes 
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and the corresponding phases are unified in the highest 
group, just as grand unified theories are unified in the 
higher groups of the symmetry breaking chain. (4) The 
unification implied by the preceding point suggests that 
the many low-lying collective states formed by system- 
atic filling of valence shells in heavy nuclei are in real- 
ity different projections in an abstract multidimensional 
space of the same state. Equivalently, the different states 
are transformed into each other by the generators of the 
symmetry. Thus, the systematics of collective modes and 
phase transitions as a function of concentration variables 
are specified by the group structure. 

It has been demonstrated that dynamical symmetries 
of the type described in Ref. Q are realized to remark- 
ably high accuracy in the spectrum and the wavefunc- 
tions of large-scale numerical calculations using the Pro- 
jected Shell Model §. Since this model is known to 
give extremely good agreement with a broad range of ex- 
perimental data (see, e.g., Ref. ||), this provides rather 
conclusive proof that these dynamical symmetries are 
strongly realized in the low-lying states of complex nuclei. 
This raises the issue of whether similar symmetries might 
be found in other complex many-body fermion systems 
such as those important in condensed matter. We may 
expect that the nuclear and condensed matter systems 
have many similarities that could make this considera- 
tion fruitful: both are composed of interacting fermions 
and both arc ultimately many-body systems that are only 
approximately describable by mean-field ideas. 

B. The Zhang SO(5) Model 

S. C. Zhang and collaborators have introduced ideas 
bearing many similarities to these into the high temper- 
ature superconductor discussion Motivated by a de- 
sire to unify AF and SC order parameters, Zhang et al 
have assembled these into a 5-dimensional vector order 
parameter, and have then constructed an 50(5) group 
that rotates the AF order parameters into the SC ones. 
This construction is based explicitly on the assumption 
of d-wave pairing in the SC state. 

In this paper we proceed differently. We start, not 
from the desire to unify two particular phases, but from 
identification of a closed algebra associated with a general 
set of fermion pairing and particle-hole operators defined 
on a periodic lattice. Nevertheless, we shall find that 
we recover Zhang's 5*0(5) symmetry as a subgroup of a 
more general U (4) symmetry if certain commutators of 
the full 5t/(4) algebra are set to zero. Thus, much of the 
extensive recent discussion of the Zhang 50(5) symmetry 
applies directly to the results of this paper. 

However, the present paper extends this discussion 
substantially: (1) The 5*0(5) subgroup is embedded in a 
larger algebra defined microscopically in the fermion de- 
grees of freedom, which implies constraints on the 50(5) 
subgroup. (2) The SU(4) highest symmetry has sub- 



groups in addition to 50(5) that may be relevant to the 
condensed matter problem in general and the cuprate su- 
perconductors in particular. (3) We frame the discussion, 
not in terms of an approximate symmetry of a Hubbard 
or t— J Hamiltonian, but in terms of an exact dynami- 
cal symmetry constructed using the Casimir invariants of 
group chains. We shall present a more detailed discus- 
sion of the relationship between the present model and 
the Zhang 50(5) model. 

II. DYNAMICAL SYMMETRY METHOD 

The dynamical symmetry method applied here corre- 
sponds schematically to the following algorithm: 

0. Assume the following conjecture: All strongly collec- 
tive modes in fermion (or boson) many-body systems can 
be put into correspondence with a closed algebra defining 
a dynamical symmetry of the sort described below. This 
is a conjecture, but there is so much evidence in support 
of it from various fields of physics that it is almost a the- 
orem: Strongly correlated motion implies a symmetry of 
the dynamics described by a Lie algebra in the second- 
quantized operators implementing that motion. 

1. Identify, within a suitable "valence space", degrees 
of freedom that one believes are physically relevant for 
the problem at hand, guided by phenomenology, theory, 
and general principles. In the present case, that reduces 
to defining a minimal set of operators that might be im- 
portant to describe superconductivity and antiferromag- 
netism on a spin lattice. 

2. Try to close a commutation algebra (of manageable 
dimension) with the second-quantized operators creating 
and annihilating the modes chosen in step 1. If necessary, 
approximate these operators, or add additional ones to 
the set if the algebra does not close naturally. In the 
present context, the simple g(k) form-factor introduced 
below is an example of simplifying things to close the 
algebra. 

3. Use standard Lie algebra theory to identify relevant 
subalgebra chains that end in algebras for conservation 
laws that one expects to be obeyed for the problem at 
hand. In the present example, we require all group chains 
to end in U(l) x SU(2), corresponding to an algebra im- 
plementing conservation of charge and spin. 

4. Construct dynamical symmetry Hamiltonians 
(Hamiltonian that are polynomials in the Casimir invari- 
ants of a group chain) for each chain. Each such group 
chain thus defines a wavefunction basis labeled by the 
eigenvalues of chain invariants (the Casimirs and the el- 
ements of the Cartan subalgebras), and a Hamiltonian 
that is diagonal in that basis (since it is constructed ex- 
plicitly from invariants). Thus, the Schoedinger equation 
is solved analytically for each chain, by construction. 

5. Calculate the physical implications of each of these 
dynamical symmetries by considering the wavefunctions, 
spectra, and transitional matrix elements of physical rel- 
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evance. This is tractable, because the eigenvalues and 
eigenvectors were obtained in step 4. Consistency of the 
symmetry requires that transition operators be related 
to group generators; otherwise transitions would mix ir- 
reducible multiplets and break the symmetry. 

6. If step 5 suggests that one is on the right track 
(meaning that a wise choice was made in step 1), one 
can write the most general Hamiltonian for the system in 
the model space, which is just a linear combination of all 
the Hamiltonians for the symmetry group chains. Since 
the Casimir operators of different group chains do not 
generally commute with each other, a Casimir invariant 
for one group chain may be a symmetry-breaking term 
for another group chain. Thus the competition between 
different dynamical symmetries and the corresponding 
phase transitions can be studied. 

7. The symmetry-limit solutions may be used as a 
starting point for more ambitious calculations that in- 
corporate symmetry breaking. Although no longer gen- 
erally analytical, such more realistic approximations may 
be solved by perturbation theory around the symmetry 
solutions (which are generally non-perturbative, so this 
is perturbation theory around a non-perturbative mini- 
mum), by numerical diagonalization of symmetry break- 
ing terms, or by coherent state or other mean-field ap- 
proximations. 

Representative application of these ideas for both 
fcrmion and boson systems may be found in nuclear 
jl],U i particle q{ , molecular , and polymer physics || . 
We also note that the general idea of symmetry having 
dynamical implications lies at the heart of local gauge 
field theories in particle physics, though the details and 
methodologies in that case differ from the ones used here. 

The only approximation in the dynamical symmetry 
approach outlined above is the space truncation. If all 
degrees of freedom are incorporated, this defines an exact 
microscopic theory. Of course, in practical calculations 
only a few carefully selected degrees of freedom can be 
included and the effect of the excluded space must be in- 
corporated by renormalized interactions in the truncated 
space. It follows that the validity of such an approach 
hinges on a wise choice of the collective degrees of free- 
dom and sufficient phcnomenological or theoretical infor- 
mation to specify the corresponding effective interactions 
of the truncated space. 



III. THE SU(4) MODEL 

Let us now introduce a mathematical formalism that 
provides a systematic implementation of the dynam- 
ical symmetry procedure for a particular physically- 
motivated choice of operators. 



A. Choice of Operators 

The success of the dynamical symmetry method de- 
pends on a wise selection of the operators that describe 
the low-energy degrees of freedom for the system. In 
the case of cuprate superconductors, we know that (un- 
like for normal superconductors) antiferromagnetism and 
superconductivity lie very near each other in the phase 
diagram. Further, data suggest that the SC phases are 
associated with Cooper pairs of spin-singlet electron holes 
having <i-wave geometry. Finally, we expect that the 
physical system must conserve both charge and spin. 
These observations suggest that at a minimum we need 
d-wave singlet pairs and operators associated with an- 
tiferromagnetism, spin operators, and charge operators 
entering the theory on an equivalent footing. Let us 
now construct a minimal closed algebra that incorporates 
these degrees of freedom |9j. 



B. The Algebra 

We begin by defining the following lattice fermion op- 
erators: 



P12 =J29( k ) c h c -ki Vxi = 



C-kiCk] 



-fc,j 9ij — (QijV 



Sij ~ c k,i Ck ,3 



(1) 



2 ^ij 



where c k ■ creates a fermion of momentum k and spin 
projection i,j = 1 or 2 = | or J,, Q = (7r,7r,7r) is an 
AF ordering vector, f2 is the lattice degeneracy, and we 
approximate a d-wave form-factor by 

g(k) = sign(cosfc x — cosk y ) = ±1 

with g(k + Q) = -g(k) and \g(k)\ = 1 (see Refs. 0,0). 
Using the usual fermion anticommutators, we deduce the 
following commutation relations among the operators of 
Eq. (Ill): 



bi2, 


Pl2] 


= —S11 — 


S22 
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ill] 


= —SikSij 
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- SkjSh 


— SjiSki 
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ill] 


= SikQij H 
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SkjQu 
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[Sij , 


Ski] 


= SjkSu - 


SuSkj 
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Qkl] 


= SjkSu - 


SuSkj 






[Sij , 


pll] 
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5 3lP\ k 
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ill] 


= s iklli + 


<Wfc 






[Qij , 


Pkl] 
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ill] 


= ^pI + 


5 3lP\k 
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Thus, this set of 16 operators is closed under commuta- 
tion and generates a Lie algebra. Detailed examination 
indicates that the algebra is associated with the group 
17(4) and has the subgroup chains 

D 30(A) x 17(1) D SU(2) S x [7(1) 
U(A) D SU(A) D SO(5) D SU(2) S x 17(1) (3) 
D SU(2) P x SU(2) S D SU(2) S x [7(1) 

that end in the subgroup SU (2) s x [7 (1) representing spin 
and charge conservation. The physical interpretation is 
aided by rewriting the generators of the U (4) algebra as 

Q+ = Qn + Q22 = X!( c fc+QT Cfe T + c k+Qi c H) 



I S12 + S21 . S12 — S21 Sn — S22 
S = [ z , -1 



2 1 2 ' 2 

Ql2 + Q21 . Ql2 — Q2I Qll ~ Q22\ 

2 ' ^ 2 ' 2 J [ ' 

# = (| (Qn - Q22), Wn + -I (Q12 + lit)) 
7?=(7?t)t rf=p\ 2 D = p 12 M = ±(n-Q) 

where Q+ generates charge density waves, S is the spin 
operator, Q is the staggered magnetization, and 7^,7? 
create and annihilate triplet d-wave pairs (see Ref. |Q), 
the operators D' and D are associated with singlet d- 
wave pairs, n is the electron number operator, and M is 
the charge operator. Properties of this group structure 
are summarized in Tables I and II, and Fig. 1. 

Notice in this context that we require exact conserva- 
tion of charge and spin for our dynamical symmetry so- 
lutions because we have not introduced approximations 
that violate these symmetries. Although it is common 
to refer to superconductivity as resulting from violation 
of particle number, this is a statement about an approx- 
imate solution. In the exact solution and in nature, par- 
ticle number is conserved |l2| . In a later section we shall 
introduce useful approximate solutions through coherent 
state methods that lead to spontaneous symmetry break- 
ing and to intrinsic states violating particle number, but 
our dynamical symmetry solutions conserve charge and 
spin exactly. 



C. The Collective Subspace 

The group SU (4) has a quadratic Casimir operator 
C su(4 ) =7r t -7r + D^D + S-S + Q-Q + M(M - 4) (5) 

The group is rank-3 and the irreducible representations 
(irreps) may be labeled by 3 weight-space quantum num- 
bers, (<ti, <T2, 03). We assume for the simplest imple- 
mentation of the model a collective d-wave pair subspace 
spanned by the following vectors: 



\S) = \n x n y n z n s ) = (n^ (7^ (ttJ (D^ |0> (6) 

This collective subspace is associated with irreps of the 
form 



(<71,Oa,<7 3 ) = (§,0,0) 



(7) 



where Q is the number of lattice sites. The correspond- 
ing expectation value of the SU (4) Casimir evaluated in 
these irreps is a constant, 



4) 



(8) 



The operator Q + defined in Eq. ([|) is the generator of 
the [7(1) factor in U{A) D U(l) x SU(A). Physically, it 
is associated with charge density wave excitations in the 
system. We note that Q + commutes with all generators 
so it annihilates the state \S), 



Q+\S) = (S\Q+\S) = 



(9) 



Thus, the collective subspace considered in isolation is 
associated with an eigenvalue Q + — 0. Physically, this 
corresponds to exclusion of charge-density wave excita- 
tions in the low-lying collective subspace of the effective 
theory || . 

The dimensionality of the full space is 2 . The di- 
mensionality of the collective subspace is much smaller, 
scaling approximately as f2 4 : 



Dim(§,0,0) = ^(§ + l)(§ + 2) 2 (§ + 3) 



(10) 



Thus for small lattices it is possible to enumerate all 
states of the collective subspace in a simple model where 
observables can be calculated analytically. 



D. SU(4) Model Hamiltonian 

The most general 2-body Hamiltonian within the d- 
wave pair space consists of a linear combination of 
(quadratic) Casimir operators C g for all subgroups g JTJ] 

g 

where Hq and H g are parameters and the Casimir oper- 
ators C g are (see Table I) 



CsO(5) 


= 7? ' ■ 7? + s ■ s - 


hM(M-3) 


CsO(4) 


= Q-Q + S-S 




C SU(2) p 


= D^D + M(M 


-1) 


CsU(2) s 


= §-S 




Cu(i) 


= M and M 2 . 





(11) 



For fixed electron number the terms in M and M 2 in Eq. 
( |ll| ) are constant. The term Hq is a quadratic function 
of particle number and may be parameterized as 



Hq = en + ^vn(n — 1), 
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where e and v are the effective single-electron energy and 
the average two-body interaction in zero-order approxi- 
mation, respectively. Thus the Hamiltonian can be writ- 
ten as 



H = H + V = en - ±vn(n - 1) 



V 



V 



-GoD^D - Givrt-vr - X Q-Q + kS-S 



(12) 
(13) 



where Go, G±, x and K are the interaction strengths of 
d-wave singlet pairing, d-wave triplet pairing, staggered 
magnetization, and spin-spin interactions, respectively. 
Since (C su ^) is a constant, by using Eq. (||) we can elim- 
inate the tv ■ tt term and after renormalizing the interac- 
tion strengths the SU(4) Hamiltonian can be expressed 



H = H' -G[ {l-p)D ] D- 
Hq = e c ff n + ±v cff n(n - 1) 



pQ-Q 



k cS S-S (14) 
(15) 



with (1 — p)G = G° ff , pG — Xeffi and k c k as the effective 
interaction strengths, and where < p < 1 for the pa- 
rameter p. Since in this paper we primarily address the 
ground state properties where 5 = 0, the last term in Eq. 
( |l4| ) will generally not enter into the later discussion. 



IV. THE DYNAMICAL SYMMETRY LIMITS 

As we have already noted, there are three subgroup 
chains of the SU(4) symmetry that conserve spin and 
charge. These define three dynamical symmetries with 
clear physical meanings. The three dynamical symme- 
try limits SU(2), 50(4), and 50(5), correspond to the 
choices p — 0, 1, and 1/2, respectively, in Eq. ([IJ). The 
Hamiltonian, eigenfunctions, energy spectrum and the 
corresponding quantum numbers of these symmetry lim- 
its are listed in Tables I and II, where we introduce a 
doping parameter x that is related to particle number 
and lattice degeneracy through 



n 



(16) 



The pairing gap A (measure of pairing order) and the 
staggered magnetization (measure of AF order) Q, 



A = G c Vz^)i/2 Q = (Q-Q) 



1/2 



(17) 



may be used to characterize the states in these symme- 
try limits. As we shall now see, each limit represents a 
different possible low-energy phase of the SU (4) system. 



A. The SO(4) Limit 

The dynamical symmetry chain 

577(4) D 50(4) x 1/(1) D SU(2) S x 17(1), 



which we shall term the SO (A) limit, corresponds to long- 
range AF order. This is the symmetry limit of Eq. (|14| ) 
when p = 1. The 50(4) subgroup is locally isomorphic to 
SU(2)p x SU(2)c, where the product group is generated 
by the linear combinations 



F=±(Q + 5) G=\{Q-S) 



(18) 



of the original 50(4) generators Q and 5. We may in- 
terpret the new generators F and G physically by noting 
that if we transform and Sy defined in Eq. (|l|) to the 
physical coordinate lattice, 



Si 



E( 

r 

= £4 



r T 
C ri C rj 



~ ^2 C ri C rj ^ 
r— even r— ode 

£ c U c n+ Yl 4 



(19) 



0(1(1 



implying that F is the generator of total spin on even sites 
and G is the generator of total spin on odd sites. Thus, 
we may interpret the 50(4) group as being generated 
by two independent spin operators: one that is the total 
spin on all sites and one that is the difference in spins 
on even and odd sites of the spatial lattice. This clearly 
is an algebraic version of the physical picture associated 
with AF long-range order. 

The 50(4) Casimir operator may be expressed as 



G so(4 ) =2(F 2 + G 2 ). 



(20) 



The 50(4) representations can be labeled by the spin-like 
quantum numbers (F — w/2, G — w/2) where w = N—fi 
with n = 0,2, ... ,N. The eigenstates are labeled by 
w and the spin 5, "0(504) = \N,w, S,m s ), and are of 
dimension (w + l) 2 . 

The ground state corresponds to u — N and 5 = 0, 
and has n/2 spin- up electrons on the even sites (F = 
N/2) and n/2 spin-down electrons on odd sites (G = 
N/2), or vice versa. Thus it has maximal staggered mag- 
netization 



|o(i -x) = h 



(21) 



and a large energy gap (associated with the correlation 
Q-Q) 



AE = 2 Xe ff(l - x)Q 



(22) 



that inhibits electronic excitation and favors magnetic in- 
sulator properties at half filling. In addition, the pairing 
gap 



A = \Gl s tl^xJT^x) 



(23) 



is small near half filling (x = 0). We conclude that these 
50(4) states are identified naturally with an AF insulat- 
ing phase of the system. 



5 



B. The SU(2) Limit 

The dynamical symmetry chain 

£17(4) D SU(2) P x SU(2) S D SU(2) S x 17(1), 

which we shall term the SU(2) limit, corresponds to 
SC order and is the p — symmetry limit of Eq. (14). 
The eigenstates are labeled by v and spin 5, tp(SU2) = 
\N, v, 5, m s ), and are of dimension (v + l)(v + 2)/2. The 
seniority-like quantum number v is the number of parti- 
cles that do not form singlet d pairs (see Table II). The 
ground state has v — 0, implying that all electrons are 
singlet-paired. In addition, there exists a large pairing 
gap 



at; = G^n 



(24) 



(see Table II) , the pairing correlation is the largest among 
the three symmetry limits, and the staggered magnetiza- 
tion vanishes in the ground state: 



= 



(25) 



Thus we propose that this state is a pair condensate as- 
sociated with a (i-wave SC phase of the cuprates. 



C. The SO(5) Limit 

The dynamical symmetry chain 

SU{4) D 50(5) D SU{2) S x 17(1), 

which we shall term the 50(5) limit, corresponds to a 
phase with the nature of a transitional or critical dy- 
namical symmetry. This symmetry limit appears when 
p = 1/2 in Eq. (Q). The SO (5) irreps are labeled by 
a quantum number t and the eigenstates may be la- 
beled by r and the spin S, ip(S05) = \N, r, S, m s ) with 
N = 0/2 — r + A, where A is the number of it pairs. The 
irreducible representation dimensionality for given N is 
(A + 1)(A + 2)/2 and the ground state has A = and 
5 = 0. 

The 50(5) dynamical symmetry has very unusual be- 
havior. Although the expectation values of A and Q for 
ground state in this symmetry limit are the same as that 
of Eq. (|2^) for the SU (2) case, there exists a huge num- 
ber of states with different values of A (the number of ir 
pairs) that can mix easily with the ground state when x 
is small because the excitation energy in this symmetry 
limit is 



AS = XG^flx 



(26) 



(see Table II). In particular, at half filling (x = 0) the 
ground state is highly degenerate with respect to A and 
mixing different numbers of it pairs in the ground state 
costs no energy. The n pairs must be responsible for the 



antiferromagnetism in this phase, since within the model 
space only ir pairs carry spin. Thus the ground state 
in this symmetry limit has large-amplitude fluctuation 
in the AF order (and SC order). This indicates that 
the 50(5) symmetry limit is associated with phases in 
which the system is extremely susceptible to fluctuations 
between AF and SC order. 



D. Energy Surfaces 

The soft nature of the 50(5) phase is seen most clearly 
if we introduce approximate solutions in terms of SU (4) 
coherent states We shall discuss the interpretation 
of 50(5) using coherent states in a separate paper |fl6| , 
but we quote one result of that study here. In Fig. 2, 
ground-state energy surfaces for various particle number 
n (or doping x) are plotted as a function of a quantity 
(3, which is related directly to the AF order parameter 
(see figure caption). Three plots are associated with the 
symmetry limits discussed above (p = 0, 1, 1/2), and one 
corresponds to a situation with a slight 50(5) symmetry 
breaking (p = 0.52). For all doping values x the energy 
minimum lies at /? = (implying that Q—0) in the SU (2) 
limit (Fig. 2a), while it lies at 







(implying that Q = n/2) for the 50(4) limit (Fig. 2b). 
In Fig. 2c, there is a broad range of doping in which 
the 50(5) energy surface is almost flat in the param- 
eter j3, implying large excursions in the AF (and SC) 
order: one can fluctuate into the other at negligible en- 
ergy cost. Notice in Fig. 2d that as doping varies the 
50(5) Hamiltonian with slight symmetry breaking inter- 
polates smoothly between AF order at half filling and 
SC order at smaller filling. Thus 50(5) acts as a kind 
of doorway between SU(2) P and 50(4) symmetries and 
this gives a precise meaning to the assertion [[| that the 
50(5) symmetry rotates between AF and SC order. 

However, the present discussion shows that the relation 
between the AF and SC phases is more complex than a 
simple 50(5) rotation because of the non-abelian nature 
of the 5/7(4) parent group of 50(5). Such dynamical 
symmetries that define a phase of the theory but that 
also interpolate between two other dynamical symmetries 
are known in nuclear structure physics where they have 
been termed critical dynamical symmetries flT^| . 

The soft 50(5) energy surface could lead to "spin- 
glass-like" phases (by which we mean phases with local 
AF or SC order but with large fluctuations in both). It 
could also lead to inhomogenous structure like stripes 
if there is a periodic spatial modulation of the system, 
since the soft nature of the energy surface implies that 
relatively small perturbations can shift an 50(5) system 
between AF and SC behavior. As noted in an earlier foot- 
note, perturbations of the symmetry by a charge density 
wave could be one source of such a spatial modulation. 
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V. PHASES AND PHASE TRANSITIONS 

In quantitative tests of the SU (4) model with parame- 
ters determined by fitting to pairing gap and pseudogap 
data, a cuprate phase diagram has been predicted. It is 
found that for the symmetry mixing parameter p close 
to but larger than 0.5 (antiferromagnetically perturbed 
50(5) symmetry) the experimental data are described 
quite well. The results of this study will be published sep- 
arately j^jnj. In this paper, we concentrate on general 
features and use the preceding discussion to construct 
the qualitative phase diagram illustrated in Fig. 3. This 
diagram, which contains essential features of the realis- 
tic phase diagram, is constructed based on the following 
arguments. 

A. Phase Diagram 

Consider the Hamiltonian ( |l4| ) and, to simplify this 
discussion, let us assume the spin term k e sS ■ S to be ne- 
glected. In that case, there are two fundamental energy 
scales in the Hamiltonian (|lj): H' and H — H' ~ Gfl 2 . 
The term H' Q originates in the single-particle energies 
e e sn and the SU(4) invariant C su ±. This term depends 
only on particle number and is isotropic: it has the same 
expectation value for any state in the S/7(4) representa- 
tion space. In contrast, the other terms H — H' (with 
characteristic energy scale Gf2 2 ) are anisotropic in the 
SU(4) space: states associated with different dynami- 
cal symmetries may have quite different expectation val- 
ues. (Thus, if we make a mean-field approximation to 
the present many-body theory — like the method of co- 
herent states — these anisotropic terms will lead to spon- 
taneously broken symmetries.) 

The term H' in the Hamiltonian ( O) may be regarded 
as the energy scale for the £7(4) D U (1) x SU(A) sym- 
metry, which describes a fermion system in which elec- 
trons are all paired but with no distinction among the 
d pairs and tt pairs. We may expect this symmetry to 
hold while the thermal energy is less than H' . When the 
system is cooled, the thermal energy eventually drops 
below the anisotropic scale Gf2 2 , the anisotropic pair- 
ing and antiferromagnetic correlations H — H' become 
important, and SU(4) breaks to its subgroups. Then dif- 
ferent low-temperature phases appear, with the favored 
phases depending on the competition between D'D and 
Q- Q interactions as a function of doping concentrations. 

From the preceding discussion, at zero temperature we 
expect the SO (4) AF phase to dominate at half filling 
if p > 0.5, because the energy surface is minimized at 
P = — x)] 1 / 2 , implying that the staggered magneti- 
zation is large: Q — n/2. On the other hand, at larger 
doping the SU (2) pairing phase is favored (the pairing en- 
ergy is optimized and the staggered magnetization min- 
imized) (see Fig. 2d). Finally, the intermediate doping 



region is described naturally by the £0(5) critical dy- 
namical symmetry that interpolates between SC and AF 
behavior with doping. Thus, ST/ (4) symmetry implies 
the schematic phase diagram of Fig. 3, independent of 
detailed calculations. 



B. Phase Transitions and Symmetry Breaking 

In the Hamiltonian (|lj) , the parameter p takes on the 
values 0, 1, and 1/2 in the three symmetry limits of the 
theory. For any other allowed value of p the system ex- 
hibits SU (4) symmetry but there is no dynamical sym- 
metry (If p ^ 0, 1, 1/2, the eigenstates of the system are 
linear superpositions of eigenstates from the three dy- 
namical symmetry chains) . In this case, phase transitions 
are driven by microscopically-determined control param- 
eters that change the expectation value of the terms of 
the Hamiltonian such that (Q-Q) is dominant in some 
circumstances while (D^-D) dominates in others. This 
permits phase transitions among the AF (SO(4) limit), 
SO(5), and SC phases (SU(2)) to be studied using a fixed 
Hamiltonian. In cuprates, the hole-doping x is an exam- 
ple of such a microscopically-determined control parame- 
ter. Thus, a Hamiltonian that possesses an approximate 
SO(5) symmetry (antiferromagnetic perturbed SO(5)) 
can have AF solutions at small hole-doping and SC so- 
lutions at large hole-doping. We shall give an explicit 
example of a phase transition driven microscopically by 
changing hole doping in the next section. 

VI. THE TRANSITION BETWEEN 
ANTIFERROMAGNETISM AND 
SUPERCONDUCTIVITY 

Above, we have used the method of generalized coher- 
ent states to give an interpretation of the SO(5) subgroup 
as a critical dynamical symmetry interpolating between 
AF and SC order and having the character of a "spin 
glass" or perhaps leading to stripe phases for a large 
range of intermediate doping parameters [ p^[ . In this 
section we address further the relationship between the 
other two phases of the theory. We show that the mi- 
croscopic symmetry incorporates a differing dependence 
on doping for SC and AF order. This implies that the 
group structure itself controls the transition between the 
superconducting SU{2) symmetry and the antiferromag- 
netically ordered SO(4) symmetry. Thus, we conclude 
that the SU (4) symmetry leads naturally to AF order at 
half-filling and to d-wave superconductivity as the system 
is hole-doped away from half-filling for a broad range of 
Hamiltonian parameters. 
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A. A Simplified Picture 

To simplify the discussion, we drop the common de- 
pendence of both phases on the spin and charge gener- 
ators and consider the competition between the £0(4) 
stabilization energy coming from a term Q ■ Q and the 
SU(2) stabilization energy coming from a term D in 
the Hamiltonian. These clearly have fundamentally dif- 
ferent behaviors with changing particle number. From 
Table II and Eq. (|l|), we may conclude that in the re- 
spective ground states, 

(xQ-Q)=xN(N + 2) 
(G D^D) = G N(Q-N + 1) 

where the pair number is N = \n. At half filling 
(N = fi/2), if x/Gq > 1, the Hamiltonian exhibits an 
effective SO (A) symmetry because the SO (A) correlation 
energy xQ ' Q dominates. As the particle number de- 
creases (hole-doping x increases), the SO(A) correlation 
energy decreases quadratically but the d-wave pairing 
D*D decreases much more slowly (essentially linearly). 
Therefore, the pairing correlation D^D will eventually 
dominate and the Hamiltonian exhibits effective SU(2) 
symmetry. 

These features imply immediately that, if x/Go > 1, 
antiferromagnetism will tend to dominate at half-filling 
but, as the system is doped away from half- filling with 
holes, eventually pairing will dominate and the system 
will become superconducting. The transition point will 
depend on the relative strengths of the D^D and Q ■ Q 
terms in the effective Hamiltonian, but the AF ground 
state at half filling and the tendency to superconductivity 
as the system is doped away from half-filling is a direct 
consequence of the group structure, independent of de- 
tailed parameter choices. 

B. Analogies in Nuclear Structure 

This competition between antiferromagnetism and su- 
perconductivity bears many strong resemblances to the 
competition in nuclear physics between spherical and de- 
formed structure for nuclei. There, it has been shown 
that the transition from spherical nuclei, which domi- 
nate the beginnings and endings of shells, to deformed 
nuclei, which often dominate the middle of shells, is 
controlled by the microscopic competition between long- 
range quadrupole-quadrupole interactions favoring defor- 
mation and short-range monopole pairing interactions fa- 
voring spherical vibrational structure. This competition 
in nuclear structure may be expressed algebraically as a 
competition between a dynamical symmetry that favors 
pairing and a dynamical symmetry that favors multipole 
(particle-hole) interactions 0. 

The essential physics of the spherical-deformed tran- 
sition in nuclear structure physics has been shown to be 



determined by the differing particle number dependence 
of the dynamical symmetries: nuclear pairing energy in- 
creases linearly with particle number but the quadrupole 
deformation energy is quadratic in particle number (that 
is, essentially the same relationship as for Eq. (E7h). 
Thus, the group structure dictates that spherical vibra- 
tional nuclei (favored by pairing energy) dominate the 
ends of shells and deformed nuclei (favored by the defor- 
mation energy) dominate the middle of shells pc| . This 
behavior is a close analog of the competition between 
antiferromagnetism dominating the half-filled lattice and 
superconductivity dominating the hole-doped lattice that 
we have discussed in this paper. 

This analogy between condensed matter and nuclear 
physics might not be accidental. The structure of heavy 
nuclei and of high-temperature superconductors both 
correspond to complex many-fermion problems in which 
strong particle-hole and pairing interactions involving d 
pairs play pivotal roles. It has been demonstrated that 
dynamical symmetries in the microscopic fermion degrees 
of freedom provide a simple but powerful unifying princi- 
ple for the nuclear structure problem |IJ] . Recent results, 
such as those presented here and in Ref. suggest 
that related symmetries may play a similar role for con- 
densed matter and that these problems from rather dif- 
ferent fields of physics may have a common dynamical 
algebraic structure. 

VII. SU(4) OR SO(5)? 

Let us now discuss more precisely the relationship be- 
tween our SU(A) model and the Zhang's SO (5) model 
j|. Although the methodologies used to obtain the two 
models are rather different, the Zhang SO(5) group is a 
subgroup of our ST/ (4) group and both model Hamiltoni- 
ans possess antifcrromagnetic perturbed SO (5) symme- 
try, implying that the two models are closely related to 
each other. Our SU(A) model and Zhang's 50(5) model 
have the same building blocks (the operator set (Q), but 
see note |22|). The essential difference is that we imple- 
ment the full quantum dynamics (commutator algebra) 
of these operators exactly, while in Ref. 0] , the dynamics 
is implemented in an approximate manner: a subset of 
10 of the operators acts as a rotation on the remaining 5 
operators {D^ , D, Q}, which are treated phenomenologi- 
cally as 5 independent components of an order-parameter 
vector (superspin n). Thus only 10 of the 15 generators 
of our SU(A) are treated fully dynamically in the Zhang 
50(5). If the full quantum dynamics (full commutator 
algebra) of the 15 operators is taken into account, the 
symmetry is SU(A), not 50(5). 

The embedding of 50(5) in our larger algebra has 
various physical consequences that do not appear if the 
50(5) subgroup is considered in isolation. We list four: 

1. As we have seen, a phase transition from an- 
tiferromagnetism to superconductivity at zero temper- 
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ature that is controlled by the doping emerges natu- 
rally and microscopically, whereas in the 50(5) model a 
symmetry-breaking term proportional to a chemical po- 
tential has to be introduced by hand. This occurs be- 
cause our 50(5) is embedded in a larger group that con- 
tains generators breaking SO(5) but preserving 577(4). 

2. The present results suggest that the 5*0(5) sub- 
group is only one of the symmetries relevant to the 
cuprate problem. It is a transitional symmetry that links 
AF to SC behavior, suggesting that it is most useful for 
the underdoped region. But the AF phases at half fill- 
ing and the optimally doped superconductors are more 
economically described by our SO (A) and SU(2) P sym- 
metries, respectively. All are unified in the SU (4) highest 
symmetry. 

3. As we discuss in a separate publication ]l8fl , the 
present SU(A) theory leads naturally to the appearance 
of pseudogap behavior (23) , which occurs above the SC 
transition temperature T c . The embedding of 50(5) in 
SU(A) is central to this behavior. 

4. The different methodology of the SU(A) dynamical 
symmetry approach suggests a different physical inter- 
pretation of the 50(5) subgroup symmetry. We suggest 
that it should be viewed as an effective symmetry that 
operates in a severely truncated space. As we shall elab- 
orate below, this implies that its microscopic validity is 
a question of the physical correctness of the matrix el- 
ements evaluated in that truncated model space, not of 
whether a particular Hamiltonian with some relevance 
for the full space (Hubbard, for example) possesses such 
a symmetry. 

Thus, both the Zhang 50(5) theory and the present 
SU{A) theory imply the existence of an approximate 
50(5) symmetry in the Hamiltonian of high-temperature 
superconductors. However, the SU(4) theory encom- 
passes a broader range of physics that may be relevant for 
cuprate superconductor. Further, as we have argued, a 
unified quantum mechanical treatment of the generators 
and order parameter vector of the Zhang SO (5) implies 
that the full quantum mechanical symmetry is SU (4) and 
not 50(5). 

VIII. DISCUSSION 

The U(A) D £7(1) x 5[/(4) symmetry represents a fully 
microscopic fermion system in which SC and AF modes 
enter on an equal footing. At this "unification" level, 
there is in a sense no distinction among these degrees 
of freedom, just as in the Standard Electroweak Theory 
of elementary particle physics the electromagnetic and 
weak interactions are unified above the intermediate vec- 
tor boson mass scale. We may expect this symmetry to 
hold while the temperature of the system is sufficiently 
high that anisotropic quantum fluctuations in the direc- 
tions associated with these collective degrees of freedom 
can be neglected (recall that the 577(4) scale is set by H' 



and the anisotropic scale by Gf2 2 ), but not so high that 
thermal fluctuations destroy the integrity of the collec- 
tive modes corresponding to the SU (4) symmetry. This 
SU(4) region then corresponds to the pseudogap regime 
& 

Generalizing language already introduced by Zhang [|4J , 
in this regime we may view the system as having con- 
densed into 577(4) pairs, which fixes the length of the 
state vectors (through the 577(4) Casimir expectation 
value) but not their direction in the state-vector space. 
Physically, this means that the system is paired, with the 
pair structure exhibiting SU(4) symmetry, but is neither 
SC nor AF because fluctuations in those directions in 
the SU(4) space are small relative to the scale set by the 
temperature. Stated in another way, the SU (4) pairs are 
of collective strength, but are not condensed into a state 
with long-range order. Stated in yet another way, neither 
the AF nor SC order parameters have finite expectation 
values in this regime, but a sum of their squares (the 
SU(4) Casimir) does. This constraint implies an inti- 
mate connection between superconductivity and antifer- 
romagnetism in the 5C/(4) model. They are, in a sense, 
different projections of the same fundamental vector in 
an abstract algebraic space. 

Compared to the Hubbard or t-J models, the dynam- 
ical symmetry approach applied here is just a different 
way to simplify a strongly-correlated electron system. In 
the Hubbard or t-J models, approximations are made to 
simplify the Hamiltonian but no specific truncation is 
assumed for the configuration space, although in prac- 
tical calculations a truncation is typically necessary. In 
contrast, we make no approximation to the Hamiltonian. 
The only approximation is the (severe) space truncation. 
The symmetry dictated Hamiltonian includes all possi- 
ble interactions in the truncated space. In principle, if 
all the degrees of freedom of the system are included, this 
approach constitutes an exact microscopic theory. The 
validity of this approach depends entirely on the validity 
of the choice of truncated space and its effective interac- 
tions, which may be tested by comparison with the data. 

Thus, we suggest that effective low-energy theories of 
the type discussed here need not have much direct rela- 
tion to the Hamiltonian or wavefunctions of the Hubbard 
or t-J models, because both the Hamiltonian and the 
wavefunctions in different model spaces could be very 
different. What is observable quantum mechanically is 
the matrix elements, not operators or wavefunctions sep- 
arately. As the SU (4) theory makes clear, we may view 
these dynamical symmetries as operating in a severely 
truncated collective subspace in which the truncation has 
been implemented primarily by symmetry considerations 
and only secondarily by energy criteria. Thus, it is pos- 
sible that the matrix elements of the S77(4) theory and 
a Hubbard or t-J model calculation might be compara- 
ble, even if the Hamiltonians and wavefunctions are sep- 
arately quite different. 

The advantage of the dynamical symmetry approach is 
its cleanness and simplicity. It is clean because the only 
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approximation is the selection of the truncated model 
space. Thus, a failure of the method is a strong signal 
that one has chosen a poor model space. It is simple be- 
cause the method supplies analytical solutions for various 
dynamical symmetry limits as a starting point. These 
symmetry-limit solutions provide an immediate handle 
on the physics and permit an initial judgement of the 
model's validity without large-scale numerical calcula- 
tions. Beyond the symmetry limits, numerical calcula- 
tions are necessary. However, because of the low dimen- 
sionality of the models spaces and the power of group the- 
ory, such numerical calculations are much simpler than 
those in, say, a Hubbard or t-J model. 

The primary limitation of the dynamical symmetry 
approach is that the interactions of the model space 
are necessarily effective interactions that will generally 
be strongly renormalized by the symmetry-constrained 
space truncation. Thus, their strengths must be sup- 
plied separately from the dynamical symmetry method 
itself, either phenomenologically (by fitting model-space 
interaction strengths to data), or by a microscopic study 
of the relationship between the effective and full-space 
Hamiltonians. 

Finally, let us counter a possible philosophical criti- 
cism of the approach taken here. Simple symmetries as a 
predictor of dynamics has found powerful application in 
fields such as particle physics or nuclear physics. How- 
ever, there is a common point of view that the possi- 
ble ground states in high-temperature superconductors 
are too complex to permit a simple model like the cur- 
rent one (which operates in a drastically truncated sub- 
space having simple wavefunctions and highly stylized 
operators) to have any validity. Although one can make 
such an argument formally, this ignores the rather ob- 
vious point that nature has managed to construct a sta- 
ble ground state having well-defined, collective properties 
that change systematically from compound to compound. 
Extensive experience in many fields of physics suggests 
that this is the signal that the phenomenon in question 
is described by a small effective subspace with renormal- 
ized interactions (that may differ substantially in form 
and strength from those of the full space) and governed 
by a symmetry structure of manageable dimensionality. 
Thus, if an approach like the one proposed here gives 
correct results for highly non-trivial phenomenology like 
the doping dependence of observable quantities, one must 
take seriously the possibility that the corresponding small 
symmetry-dictated subspace may have relevance to the 
effective behavior of real physical systems, no matter how 
complex they may appear to be superficially. 

IX. SUMMARY AND CONCLUSIONS 

In summary, an SU(4) model of High-T c supercon- 
ductivity has been proposed that contains three dynam- 
ical symmetries: A SC phase identified with the SU(2) P 



dynamical symmetry, an AF phase identified with the 
50(4) dynamical symmetry, and an S0(5) phase ex- 
tremely soft against AF and SC fluctuations over a 
substantial doping fraction that serves as a critical dy- 
namical symmetry interpolating between the other two 
phases. Realistic systems may mix these sub-symmetries 
while retaining an approximate SU(4) symmetry. Zero- 
temperature phase transitions are shown to be driven by 
the competition between the d-wave pairing and the AF 
Q • Q interactions, as controlled microscopically by the 
hole-doping concentration. This model leads naturally to 
the appearance of pseudogaps in the underdoped regime 
because it introduces multiple energy scales that permit 
pairs to form before they condense into states with long- 
range order. 

Thus, we propose that high T c behavior of the cuprates 
results from an SU(4) symmetry realized dynamically, 
and because this symmetry is microscopic its physical 
interpretation is accessible to calculation. This provides 
a solvable model incorporating many features of cuprate 
superconductors, a possible understanding of the cuprate 
phase diagram as arising from dynamically-realized sym- 
metries, and substantial insight concerning recent 5*0(5) 
theories of d-wave superconductivity. 
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FIG. 1. Dynamical symmetries associated with the U(A) 
symmetry. The generators are listed for each subgroup. 

FIG. 2. Coherent state energy surfaces for Eq. jl4]). The 
energy unit is GO, 2 /4 (see Eq. (|l4[)). The staggered magneti- 
zation Q is related to /3 by (Q) = 2Qf3 {n/2n - (3 2 ) 1/2 , where 
/3o is the value of (3 at the stable point, so /3 measures the AF 
order. Numbers on curves are the lattice occupation fractions, 
with n/Q. — 1 corresponding to half filling and < n/Q. < 1 
to finite hole doping. 



FIG. 3. Schematic phase diagram for S'[/(4) symmetry 
based on Fig. Id. The H' in Eq. (|l4|) is expressed in terms 
of hole doping x with x = 1 — n/Q; e' cS = e c gQ,/2, and Vg ff 
=v cff n 2 /2. 
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TABLE I. Properties of 5(7(4) and its subgroups.* 



Group 


Generators 


Quantum Numbers 


Casimir Operator 


Casimir Eigenvalue 


SU{4) 


Q+, S, Q, 7? f 
tt, D\ D, M 


(71 = § ((T2 = CT3 = 0) 


^■■k + D^D + S-S 
+ Q-Q + M{M -4) 


§(|+4) 


50(4) 


Q, S 


w (F — G — to/2) 


Q-Q + 5-5 


w(w + 2) 


SO(5) 


5, 7T t , TT, M 


r (w = 0) 


t^tt + 5-5 + M(M- 3) 


r(r + 3) 


SU(2) P 


D, M 


« 


D t D + M(M-l) 


i(n-D)(0-u + 2) 


SU(2) S 


5 


5 


5-5 


5(5+1) 



Valid for representations with no broken pairs. 



TABLE II. The Hamiltonian, eigenstates and spectra in three dynamical symmetry limits of the 5?7(4) model. ' 
Eg. s . is the ground state energy, AE the excitation energy, Dim the dimension of each representation for a given N, N — n/2 
the pair number, x = 1 — n/fi, , and k SQ 4 = n c ff + x e ff- 



SU(2) limit: \ip(SU2)) = \N, v, S, m s ) SO(4) limit: | ^(S04)> = \N, w, S, m s ) SO(5) limit: |j/>(S05)) = |t, JV, S, m s ) 

<C SC /(2) P > = i(n-t;)(n -t) + 2) (C SO ( 4 )) = »(ra + 2), w = N — ft <C so( 5)> = t(t + 3), r = Q/2 - JV + A 

Dim{v, N) = (v + l)(t) + 2)/2 Dim(w, N) = (w + l) 2 Dim(r, N) - (A + 1) (A + 2)/2 

H = H' - G^D^D + k cS S-S H = H' - X cff S'S + « e ff S-S H = H' a - S-S + D + K eS S-S 

(D^D) = (Csu(2) p - M(M - 1)) (S-S) = (Cso(i) - S -S) (S-S + D + D) = (C SU (4) + M - C SO (5)> 

Eg.s. =K- iG™!) 2 (l - x 2 ) E B . a . = H' — ix c ff« 2 (l ~ *f Eg.*. = H' - ixcff" 2 ! U " *) 2 

AE = vG^Q. + ft cff S(S + 1), v = w/2 A£ = MX cff (l - + k so4 5(5 + 1) AB = AG™ xQ. + K cff 5(S + 1) 

v = 0, 1,2, . . . JV; S = v,v-2,...0 = 0, 2, . . . JV; S = w, w - 1, . . . A = 0, 1, 2, . . . , N; S = A, A - 2, . . . 



Assume JV is even and 1/S7 is negligible. 
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AF Order ((3) 



Fig. 2 




Fig. 3 



